地 热 计 算 器 : 让 数值 模拟 简易 化 
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摘要 : ”地 热 作为 清洁 稳定 的 绿色 能 源 ， 在 我 国 实现 北方 清洁 供暖 全国 能 源 结构 转型 中 起 
到 重要 作用 。 可 持续 开发 地 热 资源 是 地 热 开发 利用 的 前 提 ， 而 通过 数值 模拟 的 方法 求解 数学 
模型 是 实现 地 热 资源 评价 与 提供 优化 开采 方案 以 支撑 地 热 可 持续 开发 利用 的 关键 工具 ， 得 
到 了 地 热 从 业 人 员 和 科研 工作 者 越 来 越 多 的 重视 与 关注 。 然 而 当前 并 没有 专门 针对 地 热 开 发 
利用 所 设计 的 专业 软件 ， 目 前 在 地 热 行业 使 用 较 广 泛 的 模型 软件 通常 是 相关 行业 《石油 、 水 
文 、 热 工 等 ) 软件 ， 难 以 满足 地 热 资 源 开发 利用 的 实际 需求 。 更 为 重要 的 是 ， 数 值 模拟 软件 
往往 操作 复杂 ， 使 用 极其 不 便 。 因 此 ， 我 们 开发 了 一 款 专门 针对 地 热 开 发 利用 的 中 文 软 
件 -一 “地 热 计算 器 ”， 该 软件 在 耦合 开源 数值 模拟 内 核 的 基础 上 ， 采 用 最 新 研发 的 地 热 储 
工程 优化 算法 ， 针 对 性 地 开发 了 地 热 行业 所 需 的 实际 功能 : 如 井 距 优化 、 资 源 量 评价 及 地 热 
能 开采 方案 优化 设计 等 ， 同 时 提供 友好 的 用 户 交 互 界面 ， 实 现 软件 的 “简易 化 ”或 者 “一 
键 操 作 ” 功 能 ， 证 具备 不 同 专业 背景 的 地 热 从 业 人 员 均 能 够 在 短 时 间 内 掌握 其 操作 ， 一 方 
面 促进 地 热 行 业 数值 模拟 技术 的 应 用 ， 另 一 方面 也 为 地 热 资 源 的 可 持续 开发 利用 提供 服务 。 
关键 词 : 地 热 计 算 器 ;数值 模拟 ; 地 热能 ， 可 持续 开发 利用 
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Abstract: 
Geothermal energy plays a key role in helping promote the clean space heating 


in North China and National energy structure optimization. The sustainable 
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geothermal energy utilization needs the numerical simulation to solve the 
issues on geothermal resources assessment, strategy optimization and etc. ， 
which attracts more and more attention. However, there are no numerical 
simulation software specializing on geothermal energy exploitation and 
utilization. As a result, the geothermal managers have to use the numerical 
software in other fields (0il, Hydrology and Thermal Engineering etc.). What” s 
more, the existing software is always hard to operate. Here in this paper, we 
develop a new software named “Geothermal kits” , which is developed by 
coupling the geothermal optimal algorithm and open-source numerical simulation 
software. The ‘Geothemral kits’” software has very friend interface, which 
could be done by one click. Therefore, it could be easily mastered by the 
geothermal scientists, managers and engineers with different professional 
background. The ‘Geothemral kits” will help solving the problems related to 
numerical simulation and serve for the sustainable geothermal energy 
utilization. 
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数学 模型 是 可 持续 开发 地 热 资 源 的 必 备 工具 。 从 地 热 资源 的 勘查 评价 到 地 热 
资源 开发 利用 的 优化 方案 制定 ， 都 离 不 开 数 学 模型 。 数 学 模型 的 运算 主要 包括 解 
析 法 和 数值 法 两 种 。 一 般 而 言 解析 法 简单 、 运 算 快 并 且 能 够 获得 精确 的 解 。 但 在 多 
数 与 地 热 相 关 的 问题 中 ， 要 考虑 地 质 和 地 热 地 质 条 件 的 复杂 性 ， 解 析 解 往往 无 
符 。 因 此， 数值 法 即 通 过 数值 模拟 进行 运算 ， 成 为 我 们 在 地 热能 开发 利用 过 程 中 
最 常用 的 解法 。 事 实 上 ， 在 2020 年 5 月 1 日 起 实施 的 地 质 矿产 行业 标准 《DZ/T 
0331-2020: 地 热 资源 评价 方法 及 估算 规程 以 及 尚 在 征求 意见 中 的 国家 标准 《地 
源 热泵 系统 工程 技术 规范 征求 意见 稿 ) 均 对 资源 勘查 程度 较 高 的 地 热 田 ， 推 荐 
采用 数值 模拟 的 方法 评价 资源 量 以 及 制定 优化 的 开采 方案 。 

数值 模拟 是 地 热学 乃至 地 热 资 源 开发 利用 研究 中 的 重要 工具 ,在 地 球 动力 学 
研究 中 ， 采 用 数值 模拟 热 演 化 规律 是 其 重要 研究 内 容 之 一 趾 ， 针 对 中 深层 地 热 
资源 的 开发 ， 可 以 采用 数值 模拟 方法 给 出 采 热 上 限 钼 ， 在 干 热 岩 开发 乃至 超 临 
界 地 热流 体 的 开发 ， 数 值 模拟 均 可 完成 室内 实验 无 法 模拟 的 物理 过 程 54。 以 

“Geothermal numerical modelling (地热 数值 模拟 ) ”为 检索 关键 词 ， 在 Web 
of Science 网 站 进行 检索 ， 可 以 看 到 过 去 30 年 来 ， 无 论 是 发 文 量 还 是 引文 量 都 
在 快速 增加 ， 显 示 地 热 数 值 模拟 的 被 关注 度 以 及 影响 力 在 逐年 提高 (图 1) 。 
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图 1 以 “地 热 数值 模拟 ”为 检索 关键 词 在 Web of Science 网 站 
(http://apps.webofknowledge.com) 检索 到 的 论文 量 和 相应 被 引 频次 


目前 国际 上 常用 的 地 热 数 值 模拟 软件 包括 Comsol、Feflow、0penGeoSys 和 
Tough2 等 。 当 然 ， 以 上 软件 更 类 似 于 “计算 器 ”的 功能 ， 它 们 不 仅 可 以 用 在 地 
热能 的 开发 利用 上 ， 在 水 污染 、 核 废物 处 置 和 二 氧化 碳 封存 等 诸多 领域 均 可 采用 
近 些 年 ， 国 内 科研 界 广泛 采用 上 述 软件 ， 提 出 了 大 量 用 于 地 热 开 发 比如 资源 评 
价 、 井 距 评价 和 优化 开采 等 的 研究 算法 。 但 令 人 遗憾 的 是 ， 在 国内 的 地 热 产 业界 
数值 模拟 以 及 基于 数值 模拟 技术 的 上 述 算法 的 应 用 十 分 有 限 。 其 原因 一 方面 是 上 
述 软件 以 英语 为 主 、 操 作 较 为 繁琐 ， 短 时 间 内 难以 掌握 ， 另 一 方面 是 软件 收费 多 
较为 昂贵 ， 应 用 成 本 较 高 。 

针对 以 上 问题 ， 我 们 开发 了 一 球 中 文 软件 ， 将 之 命名 为 “地 热 计 算 器 ”， 
以 实现 基于 数值 模拟 软件 的 计算 功能 ， 并 将 最 新 研发 的 新 算法 耦合 于 软件 之 中 ， 
力争 实现 软件 的 “简易 化 ”或 者 “一 键 操 作 ” 功 能 ， 让 有 具备 不 同 专业 背景 的 从 
业 人 员 均 能 够 在 短 时 间 内 掌握 其 操作 ， 以 促进 地 热 行 业 数值 模拟 技术 的 应 用 ， 
为 地 热 资 源 的 可 持续 开发 利用 提供 服务 。 


2 软件 介绍 


地 热 计算 器 是 由 Python 语言 编写 ， 它 的 主要 功能 是 通过 调用 数值 模拟 软件 ， 
实现 为 地 热 资源 可 持续 开发 服务 。 地 热 计算 器 的 初步 构架 包括 井 距 优化 、 资 源 评 
价 、 中 深层 地 埋 管 换 热 量 评估 、 地 热 开发 优化 方案 设计 等 模块 。 当 前 1. 0 版 本 仅 能 
够 实现 井 距 优化 功能 ， 因 此 本 文 着 重 介绍 井 距 优化 模块 。 

在 地 热能 的 开发 利用 过 程 中 ， 需 通过 回 灌 以 实现 可 持续 利用 ， 但 回 灌 井 和 
生产 井 相 近 ， 易 发 生 热 突破 ， 二 者 相距 太 远 ， 又 难以 维持 热 储 压力 。 为 实现 井 距 


的 优化 设计 ，Kong et al. (2017) 提 出 了 面向 经 济 成 本 的 优化 算法 外， 该 算法 
主要 分 为 两 步 : 第 一 步 是 基于 热 田 的 地 热 地 质 概 念 模型 ， 建 立 耦 合 地 下 水 流 和 
热传导 两 个 过 程 的 数值 模型 ， 输 出 不 同 采 滩 情 景 下 的 温 压 响应 ;， 第 二 步 是 根据 
温 压 响应 结果 计算 由 此 而 带 来 的 经 济 损失 ， 从 而 以 经 济 最 优 为 目标 函数 ， 给 出 
最 优 采 灌 布局 名 。 

“地 热 计算 器 ”软件 井 距 计算 的 核心 算法 基于 Kong et al. (2017) 四， 计 
算 过 程 中 调用 了 德国 辫 姆 霍 兹 环境 研究 中 心 0laf Kolditz 团队 开发 的 开源 数值 
模拟 软件 0OpenGeoSys""。 本 文 着 重 介 绍 “ 地 热 计 算 器 ”的 编写 原理 ， 有 关 井 距 
优化 算法 和 0penGeoSys 的 详情 请 分 别 参 加 Kong et al. (2017) 中 和 Kolditz 
et al.，201206] 。 

图 2 为 “地 热 计算 器 ”的 计算 流程 : 在 打开 软件 后 ， 需 要 用 户 首先 输入 相 
应 的 参数 ， 包 括 储 层 与 开发 利用 的 相关 参数 等 ， 可 采用 用 户 自 定义 或 默认 人 缺 省 
值 两 种 模式 ， 程 序 会 自动 判别 用 户 输入 的 参数 是 否 合理 ， 然 后 程序 将 合理 的 参 
数 传递 给 0penGeoSys 计算 内 核 ， 自 动 构建 输入 文件 ， 并 调用 0penGeoSys 可 执 
行文 件 执行 计算 ; 待 计 算 完成 后 程序 将 自动 执行 模型 后 处 理 ， 分 别 输出 不 同 井 
距 下 的 地 温 场 、 压 力 场 的 分 布 ， 并 基于 此 计算 不 同 井 距 下 的 开采 成 本 ， 最 后 综合 
相关 信息 ， 程 序 给 出 最 优 井 距 的 推荐 值 。 


计算 每 年 分 别 由 温度 和 压力 降低 引起 的 开采 成 本 增加 
+ 


取 成 本 最 低 情况 下 的 井 距 为 最 优 井 距 ， 并 绘图 
} 


图 2 地 热 计算 器 软件 计算 流程 图 


3 ”案例 介绍 


本 文 以 我 国 北方 某 地 区 碳酸 盐 岩 地 热 田 的 地 热 资 源 开 发 为 例 ， 介 绍 地 热 计 
算 器 的 计算 结果 。 该 含水 层 厚 300m， 储 层 类 型 为 碳酸 盐 岩 ， 孔 阶 度 和 渗透 率 分 
别 为 0.3 和 1X102 亚 ;， 该 地 区 每 年 供暖 4 个 月 ， 其 余 时 间 停 采 ， 生 产 井 流量 
150 mh， 生产 井 初 始 水 温 为 80'C ，100% 回 洪 ， 回 灌水 温 为 35'C。 给 定 热 田 寿 
命 50 年 ， 构 建 该 地 区 的 数值 模型 ， 所 需 参 数 如 表 1. 


表 1: 地热 计算 器 输入 参数 表 

基 了 筷 渗 从 生 9 开采 ”岩石 热 。” 岩石 比 热 ” 岩石 ”开采 回 
础 际 透 水 产 滩 着 前 导 率 容 密度 流量 灌 
参 度 率 层 井 水 温 水 位 流 
数 厚 出 量 
与 度 上 书 
取 水 
值 温 

0.3 1x 300 80 35 0 2.8 950 2500 150 150 

10B 
m m C CC m W/m* Kk) J/(kg*k) Kg/m m/h mh 

高 地 区 贮 热 ， 热 价 电价 贴现 率 
级 热 域 水 而 
参 田 地 率 
数 寿 下 效 
与 命 水 率 
取 流 
值 速 

50 5Xx Dx 0.7 88.9 0.58 0. 046 

10°5 10° 
年 m/s ml! ¥/GJ] ¥/kWh 


考虑 到 地 质 学 模型 和 经 济 学 模型 中 的 部 分 参数 不 常见 ，“ 地 热 计算 器 ” 软 
件 做 了 “高 级 设置 ”功能 ， 在 用 户 对 高 级 参数 概念 不 清晰 的 背景 下 ， 可 不 予 调 
整 这 些 参数 。 

基于 以 上 参数 ， 通 过 单 击 “ 计 算 并 绘图 ”按钮 ， 便 可 实现 一 键 计 算 功 能 。 地 
热 计算 器 的 计算 结果 以 图 和 动画 的 形式 呈现 ， 分 别 如 图 3- 图 6 所 示 。 
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3 生产 井 水 位 同 井 距 的 关系 


图 3 给 出 了 热 田 开采 50 年 后 不 同 井 距 下 生产 井 的 水 位 变化 。 从 中 可 以 看 出 ， 
在 井 距 为 100m 时 ， 水 位 与 原始 水 位 相 比 下 降 显 著 ， 但 是 在 200m 之后， 水 位 下 
降 量 不 再 随 井 距 发 生变 化 ， 说 明 此 时 ， 将 回 灌 井 设置 更 远 ， 对 水 位 已 无 显著 意 
义 。 
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4 生产 井 温 度 与 井 距 的 关系 
图 4 给 出 了 50 年 后 生产 井 水 温 同 井 距 的 关系 。 可 以 看 到 ， 在 井 距 小 于 300m 


时 ， 生 产 井 水 温 下 降 最 著 ， 到 井 距 为 400m 时 ， 才 避免 了 热 突破 的 发 生 。 而 后 ， 
随 着 井 距 的 加 大 ， 生 产 井 水 温 基 本 不 再 下 降 。 
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图 5 经 济 成 本 与 井 距 的 关系 
5 给 出 了 50 年 内 由 于 生产 井 水 位 下 降 和 水 温 下 降 而 引起 的 成 本 损失 随 井 
距 的 变化 规律 。 可 以 看 到 当 井 距 大 于 400m 时 ， 经 济 成 本 已 不 再 发 生 显著 变化 。 
“地 热 计 算 器 ”软件 通过 自动 求 取 不 同 井 距 下 的 经 济 成 本 ， 获 得 最 优 井 距 为 
400m。 对 比 图 3 和 图 5 可 以 发 现 ， 经 济 成 本 和 生产 井 水 温 的 变化 显著 相关 ， 这 也 
山 显 了 热 田 生产 应 尽 可 能 降低 生产 井 水 温 下 降 的 重要 性 。 


chinaXiv:202007.00048vVv2 


温度 场 分 布 图 (第 32 年 ， 景 优 井 距 为 400m) 
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水 流 场 分 布 图 (第 32 年 ， 最 优 并 距 为 400m) 


6 热 储 层 水 流 场 与 温度 场 实 时 分 布 
图 6 在 “地 热 计 算 器 ”软件 界面 内 以 “动画 ”形式 展示 ， 它 形象 的 展示 了 
热 储 层 温 度 和 压力 随时 间 发 生 的 变化 。 由 于 模型 设置 的 步 长 为 1 年 1 个 周期 ， 所 


以 该 动画 的 展示 频率 同样 以 年 为 单位 进行 展示 ， 辅 助 地 热 资源 开发 的 决策 者 进 
行 热 田 生产 管理 。 


4 结论 与 前 瞻 


“地 热 计算 器 ”软件 基于 地 热 可 持续 开发 利用 的 最 新 算法 ， 通 过 调用 数值 模 
拟 计算 内 核 ， 一 键 实现 与 地 热 资源 开发 利用 的 相关 功能 ， 以 使 得 数值 模拟 在 地 
热 行业 的 应 用 更 加 “简易 化 ”/“ 傻 瓜 化 ”。 软 件 输出 的 结果 以 图 片 和 动画 的 形 
式 展示 ， 便 于 从 业 人 员 直 接 拿 来 使 用 。 尽管 当前 软件 版 本 (V1.0) 仪 包含 有 井 距 
优化 功能 ， 在 未 来 的 升级 版 本 中 会 逐步 加 入 以 下 功能 : 

(1) 资源 评价 。 该 模块 含 基础 资源 量 评价 和 动态 资源 评价 。 基 础 资源 评价 将 在 
三 维 地 质 模型 的 基础 之 上 完成 ， 提 高 评价 精度 ， 以 打破 简单 的 体积 法 计算 ; 动 
态 资源 量 评价 将 考虑 不 同 的 生产 模式 ， 直 接 服 务 于 生产 。 

《2) 中 深层 地 源 热泵 换 热量 评估 。 该 模块 将 计算 给 定 地 热 地 质 条 件 下 的 换 热 
量 ， 并 分 析 不 同 换 热量 下 岩 土 的 温 压 响应 。 

(3) 地 热 优 化 开采 方案 设计 。 该 模块 将 计算 给 定 地 热 条 件 下 ， 最 优 的 开采 模 
式 ， 包 括 开采 量 、 回 灌 量 、 回 灌水 温 等 一 系列 参数 ， 并 将 完成 群 井 采 灌 模 式 下 的 
综合 优化 模式 设计 。 

(4) 示 踪 试验 结果 分 析 。 该 模块 将 针对 示 踪 试验 的 监测 结果 ， 一 键 给 出 渗透 
系数 、 地 下 水 流速 和 热 突破 预测 等 关键 信息 。 

(5) 浅 层 地 热能 模块 。 该 模块 将 会 围绕 浅 层 地 热能 开发 的 一 系列 需求 ， 耦 合 
地 源 热泵 模块 ， 实 现 简易 化 操作 。 
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